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Abstract 

We present the first non-trivial online algorithms for the non-uniform, multicommodity buy-at-bulk 
(MC-BB) network design problem. Our competitive ratios qualitatively match the best known approxi¬ 
mation factors for the corresponding offline problems. In particular, we show 

• A polynomial time online algorithm with a poly-logarithmic competitive ratio for the MC-BB 
problem in undirected edge-weighted graphs. 

• A quasi-polynomial time online algorithm with a poly-logarithmic competitive ratio for the MC- 
BB problem in undirected node-weighted graphs. 

• For any fixed e > 0, a polynomial time online algorithm with a competitive ratio of 
(where k is the number of demands, and 0(.) hides polylog factors) for MC-BB in directed graphs. 

• Algorithms with matching competitive ratios for the prize-collecting variants of all the above prob¬ 
lems. 

Prior to our work, a logarithmic competitive ratio was known for undirected, edge-weighted graphs 
only for the special case of uniform costs (Awerbuch and Azar, FOCS 1997), and a polylogarithmic 
competitive ratio was known for the edge-weighted single-sink problem (Meyerson, SPAA 2004). To 
the best of our knowledge, no previous online algorithm was known, even for uniform costs, in the 
node-weighted and directed settings. 

Our main engine for the results above is an online reduction theorem of MC-BB problems to their 
single-sink (SS-BB) counterparts. We use the concept of junction-tree solutions (Chekuri et al, FOCS 
2006) that play an important role in solving the offline versions of the problem via a greedy subroutine 
- an inherently offline procedure. Our main technical contribution is in designing an online algorithm 
using only the existence of good junction-trees to reduce an MC-BB instance to multiple SS-BB sub¬ 
instances. Along the way, we also give the first non-trivial online node-weighted/directed single-sink 
buy-at-bulk algorithms. In addition to the new results, our generic reduction also yields new proofs of 
recent results for the online node-weighted Steiner forest and online group Steiner forest problems. 
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1 Introduction 


In a typical network design problem, one has to find a minimum cost (sub) network satisfying various 
connectivity and routing requirements. These are fundamental problems in combinatorial optimization, 
operations research, and computer science. To model economies of scale in network design, Salman et 
al. [32] proposed the buy-at-bulk framework, which has been studied extensively over the last two decades 
(e.g., [6, 20, 33, 22, 29, 28, 15, 14]). In this framework, each network element is associated with a sub¬ 
additive function representing the cost for a given utilization. Given a set of connectivity demands com¬ 
prising k source-sink pairs, the goal is to route integral flows from the sources to the corresponding sinks 
concurrently to minimize the total cost of the routing. 

An important application of the problem is capacity planning in telecommunication networks or in the 
Internet. As observed by Awerbuch and Azar [6], this application is inherently “online” in that terminal- 
pairs arrive over time and need to be served without knowledge of future pairs. The authors of [6] give 
a logarithmic-competitive online algorithm for the uniform case where every edge is associated with the 
same cost function. However, uniformity is not always a feasible assumption, especially in heterogeneous, 
dynamic networks like the Internet. Indeed recent research (e.g., [29, 28, 15]) has focused on the non- 
uniform setting with a different sub-additive function for every network element. In this non-uniform setting, 
Meyerson [28] gives a polylogarithmic-competitive algorithm for the special case when all terminal-pairs 
share the same sink. To the best of our knowledge, no non-trivial online algorithm is known for the general 
multicommodity setting, which is the focus of our paper. 

We consider, in increasing order of generality, undirected edge-weighted graphs, undirected node- 
weighted graphs, and directed edge-weighted graphs.^ It is also convenient to classify the problems that 
we study into the single-sink version where all the terminal-pairs share a common sink, and the general 
multicommodity version where the sinks in the terminal-pairs may be distinct. For notational convenience, 
we use the following shorthand forms for our problems: X-Y-BB where X = SS or MC (single-sink and 
multicommodity, respectively) and Y = E or N or D (undirected edge-weighted, undirected node-weighted, 
and the general directed case, respectively). 

1.1 Our Contributions 

We obtain the following new results (unless otherwise noted, our algorithms run in polynomial time): 

• A poly-logarithmic competitive online algorithm for the MC-E-BB problem. 

• A poly-logarithmic competitive online algorithm for MC-N-BB and SS-N-BB that runs in quasi¬ 
polynomial time. 

• An 0(/c 2 +^)-competitive online algorithm for MC-D-BB for any constant e > 0 with running time 

where 0{.) hides polylogarithmic factors. Eor SS-D-BB, the ratio improves to 0{k^), trans¬ 
lating to a polylogarithmic competitive ratio in quasi-polynomial time. 

• Online algorithms for prize-collecting versions of all the above problems with the same competitive 
ratio. 

Up to exponents in the logarithm, our online algorithms match the best known offline approximation al¬ 
gorithms (Chekuri et al. [15] for MC-E/N-BB and Antonakopoulos [5] for MC-D-BB); for MC-N-BB, 

*In undirected graphs, node costs can simulate edge costs; in directed graphs they are equivalent. 
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however, a polynomial time, polylogarithmic approximation is known [14], whereas our algorithm runs in 
quasi-polynomial time. Furthermore, a logarithmic lower bound, even for SS-E-BB, follows from the lower 
bound for the online Steiner tree problem [26], and a polylogarithmic lower bound for online SS-N-BB 
follows from a matching one for set cover [3]. 

From a technical perspective, we derive all the multicommodity results using a generic online reduction 
theorem that reduces a multicommodity instance to several single-sink instances, for which we either use 
existing online algorithms or give new online algorithms. Informally, one can view this as the “online 
analog” of the junction-tree approach pioneered by Chekuri et al. [15] for offline multicommodity network 
design. We discuss this approach in the next subsection. 

1.2 An Online Reduction to Single Sink Instances 

Multicommodity network design problems, both online and offline, are typically more challenging than their 
single-sink counterparts, and have historically^ required new ideas every time depending on the specific 
problem at hand. The situation is no different for buy-at-bulk, both for uniform and non-uniform costs. 

In the offline buy-at-bulk setting, this shortcoming is addressed by Chekuri et al. [15] (expanded to 
other problems by [14, 13, 5]), who introduce a generic combinatorial framework for mapping a single 
instance of a multicommodity problem to multiple instances of the corresponding single-sink problem. At 
the heart of this scheme is the following observation that holds for many multicommodity problems such as 
(edge/node) Steiner forest, directed Steiner network, buy-at-bulk, and set connectivity: there exists a near- 
optimal^ junction-tree solution for the multicommodity problem that decomposes into solutions to multiple 
single-sink problems where each single-sink problem connects some subset of the original terminal-pairs to 
a particular root. 

The problem now reduces to finding good junction-trees to cover all the terminal-pairs. The offline tech¬ 
niques [15, 13, 5] tackle this using a greedy algorithm for finding the single-sink solutions; more precisely, 
in each step they find the best density (cost per terminal-pair) solution that routes a subset of terminal- 
pairs via a single sink. A set cover style analysis then bounds the loss for repeating this procedure until 
all terminal-pairs are covered. However, as the reader may have already noticed, the greedy optimization 
approach is inherently offline, as finding the best-density solution requires knowledge of all terminal-pairs 
upfront. Our main technical contribution in this work is an online version of the junction-tree framework. 
Indeed, we show how to reduce any multicommodity buy-at-bulk instance to a collection of single-sink 
instances online. 

(Informal Theorem) If the junction-tree approximation factor of the MC-BB problem is a, the integrality 
gap of a natural LP relaxation of the SS-BB problem is /3, and there is a y-competitive online algorithm for 
the SS-BB problem, then there is an 0{al3y ■ polylog(n))-competitive algorithm for the MC-BB problem. 

To prove the above theorem, we first write a composite-LP, which has (a) an outer-LP comprising assign¬ 
ment variables that fractionally assign terminal-pairs to roots, and (b) many inner-LPs which correspond to 
the natural LP relaxations for the SS-BB problem for each root and the terminal-pairs fractionally assigned 
to it by the outer-LP. We then apply the framework of online primal-dual algorithms (see [10] for instance) 
to solve the composite-LP online. However, there are two main challenges we need to surmount. 

• First, the existing framework has been mostly applied to purely covering/packing LPs^ and our inner- 

^For instance, compare [29] and [12] for the SS-E-BB and MC-E-BB problem, and compare Naor et at. [31] and Hajiaghayi et 
al. [24] for the online node-weighted Steiner tree and Steiner forest problem. 

^We call the quality of such a solution the junction-tree approximation factor; e.g., it is O(logn) for MC-E-BB and MC-N- 
BB [15] 

"'Our current understanding of mixed packing-covering is limited [7] and does not capture the problem we want to solve. 
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LPs have both kinds of constraints, and moreover, there is an outer-LP encapsulating them. We show 
nonetheless that it can be extended to solving our LP fractionally up to polylogarithmic factors. Indeed, 
we use the specific flow-structure of the inner-LP, and each step of our algorithm solves many (auxiliary) 
min-cost max-flow problems. 

• The second difficulty is in rounding this fractional solution online. This is a hard problem, and cur¬ 
rently we do not know how to do so even for basic network design problems such as the Steiner tree problem. 
To circumvent this, we show that it suffices to only partially round the LP. More precisely, we round the 
LP solution so that only the outer-LP (assignment variables) become integral, and the inner-LPs remain 
fractional. This gives us an integral assignment of the terminal-pairs to different single-sink instances, with 
bounded total fractional cost. Now, from the bounded integrality gap of the inner-LPs, we know that there 
exist good single-sink solutions for our assignment of terminal-pairs to roots, even though we cannot find 
them online^. Using this knowledge, our final step is to run online single-sink algorithms for each root, 
and send the terminal-pairs to the root as determined by the outer-LP assignment. Figure 1 summarizes our 
overall approach. 


initialize multiple online algorithms for the single-sink problem (one for each vertex as root), 
when {si,ti) arrives 

update the fractional solution of the composite LP to satisfy the new request 

round the composite LP to get an integral solution to the outer LP, which gives us an assignment 
of {si,ti) to some root r 

send both Si and ti to the instance of the single-sink online algorithm with root r 


Figure 1: Online Framework for Multicommodity Network Design Problems 

The results mentioned in Section 1.1 follow by bounding a, /?, 7 for the corresponding problems. For 
MC-E-BB, all of these are known to be 0(polylog(n)) ([15, 29, 28] respectively). For MC-N-BB, it is 
known both a, /3 are bounded by 0(polylog(n)) [14], and we bound 7 in Section 5 by giving the first online 
algorithms for SS-N-BB. For MC-D-BB, we need some additional work. In this case we cannot directly 
bound /3, since the integrality gap of the natural LP relaxation is not known to be bounded. Nevertheless, 
in Section 4, we show that it suffices to work only with more structured instances for which we can bound 
the integrality gap. 

Finally, we illustrate the generality of our reduction theorem by noting that, when combined with existing 
bounds on a, /?, and 7 , it immediately implies (up to polylogarithmic factors) some recent results in online 
network design, such as online node-weighted Steiner forest [24], and online edge-weighted group Steiner 
forest [31] - two problems for which specialized techniques were needed, even though their single-sink 
counterparts were known earlier. 

^The difficulty comes from the fact that the online solution we maintain must be monotonic, i.e., the decisions are irrevocable. 
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1.3 Related Work 

Buy-at-bulk network design problems have reeeived eonsiderable attention over the last two deeades, both 
in the offline and online settings. For the uniform eost model, Awerbueh and Azar [6] give an O(logn)- 
approximation for MC-E-BB, while 0(l)-approximations are known [20, 33, 22] for SS-E-BB. We also 
note that 0(l)-approximations have been obtained in speeial eases for the multieommodity problem, sueh 
as in the rent-or-buy setting [21]. Meyerson et al. [29] give an 0(logA:) approximation for the general 
SS-E-BB, and the first non-trivial algorithm for MC-E-BB is an exp(0(-v/log n log log n))-approximation 
due to Charikar and Karagiozova [12]. This was improved to a poly-logarithmie faetor by Chekuri et al. 
[15] who also solve MC-N-BB [14] with similar guarantees. Eor direeted graphs, our knowledge is mueh 
sparser. Even for speeial eases like direeted Steiner tree and forest, the best polytime approximation faetors 
known are 0{k^) [34, 11] and min(0(\/fe, [13, 18, 8] respeetively, and these ideas were extended 

to MC-D-BB by Antonakopoulos [5]. On the hardness side, Andrews [4] shows that even the MC-E-BB 
problem is 0( log^/^“^ n)-hard, while MC-D-BB (in faet direeted Steiner forest) is known to be label-eover 
hard [17]. 

The online Steiner tree problem (a speeial ease of online SS-E-BB) was first studied by Imase and 
Waxman [26] who give an 0(log A:)-eompetitive algorithm. Berman and Coulston [9] give an 0(log/c)- 
eompetitive algorithm for online Steiner forest, and both these results are tight, i.e., there is an 0(log k) lower 
bound. As mentioned earlier, Awerbueh and Azar’s algorithm [6] ean be seen as an 0(logn)-eompetitive 
online algorithm for the uniform-cost MC-E-BB. Eor non-uniform buy-at-bulk, the only online algorithm 
that we are aware of is Meyerson’s [28] polylog-eompetitive algorithm for the single-sink problem. Eor 
online node-weighted network design, developments are mueh more reeent. A polylogarithmie approxima¬ 
tion for the node-weighted Steiner tree problem was first given by Naor et al. [31] and later extended to the 
online node-weighted Steiner forest problem [24] and prize-eolleeting versions [23]. These algorithms, like 
ours in this paper, utilize the online adaptation of the primal-dual and EP rounding sehemas pioneered by 
the work of Alon et al. [3] for the online set eover problem (see also [2] for its adaptation to network design 
problems). We also note that in the node-weighted setting, the online lower bound ean be strengthened to 
D(log n log k) using online set eover lower bounds [3, 27]. 

2 Preliminaries and Results 

We now formally state the problem, set up notation that we use throughout the paper, and state our main 
theorems. 

2.1 Our Problems 

Buy-at-bulk Network Design. In the most general setting of the MC-D-BB problem, an instanee X = 
{G,X) eonsists of a direeted graph G = {V,E) and a eolleetion X of terminal-pairs {si,ti) € V x V', 
eaeh sueh Si and U is ealled a terminal. Eaeh pair also has a positive integer demand di, whieh we 

assume to be 1 for elarity in presentation.^ Additionally, eaeh edge e G 77 is assoeiated with a monotone, 
sub-additive^ eost funetion /e ; M>o —)• lR>o- A feasible solution to the problem is a eolleetion of paths 
{Pi,..., Pk} where Pi is a direeted path from Si to f earrying load di. Given a solution {Pi,..., P^}, we 

®We can handle non-uniform demands by incurring an additional 0(log D) factor in the competitive ratio and the running time 
(where D is the maximum demand) by having 0{logD) “unit-demand” instances, where the instance deals with demands 
between 2®“^ and 2*. 

’That is, fe{x) > fe{y) whenever x > y, and fe{x y) < fe{x) -f fe{y) 
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let load(e) = Ei-.e&Pi di denote the total load on edge e. The goal is to find a feasible solution minimizing 
the objective Objgg := Yle&E /e(load(e)). In the online problem, the offline input consists of the graph G 
and the cost functions /g. The pairs (sj, ti) arrive online in an unknown, possibly adversarial, order. When a 
pair [si, ti) arrives, the algorithm must select the path Pi that connects them, and this decision is irrevocable. 

Reduction to the Two-metric Problem. Following previous work, throughout this paper we consider an 
equivalent problem (up to constant factors) known as two-metric network design. In this problem, instead 
of functions /e(.) on the edges, we are given two parameters Cg and £g on each edge. One can think of 
Cg as a fixed buying cost, or just cost, of edge e, and as a per-unit flow cost, or length, of edge e. The 
feasible solution space is the same as for the buy-at-bulk problem, and the goal is to minimize the objective 
0 bj2M := ZleeUi P ^eeP; following lemma is well known (see e.g., [15]). 

Lemma 1. Given an instance of the buy-at-bulk problem, for any e > 0 one can find an instance of the 
two-metric network design problem such that, for any feasible solution, Obj 2 i\/| < Objgg < (2 + e) 0 bj 2 M- 

Remark. In light of the above lemma, henceforth we abuse notation and let the buy-at-bulk problem mean 
the two-metric network design problem. 

2.2 Our Tools 

Junction-tree solutions. Given an instance Z = {G, X) of the buy-at-bulk problem, we consider junction- 
tree^ solutions, a specific kind of solution to the problem introduced by [15]. In such solutions, the collection 
of pairs are partitioned into groups and each group is indexed by a root vertex r G fo. For all terminal pairs 
{si,ti) in a group indexed by r, the path Pi from Si to ti contains the root vertex r (see Figure 2). 



Figure 2: A group of terminal pairs routed via a junction-vertex in an undirected graph. 

Formally, consider an instance X = {G, X) of the buy-at-bulk problem, and let Opt denote the objective 
value of the optimum solution. Given a partition Ft := (tt^i , • • •, itr ^) of terminal pairs indexed by q different 
root vertices, a junction-tree solution is one that uses single-sink solutions to connect the original terminal- 
pairs. Indeed, for each part iTr indexed by root r, consider the optimal solutions to the single-sink problem 
on graph G with demands {(sj,r) : {si,ti) G vr^} and the single-source problem^ with pairs {(r, tj) : 
{si,ti) G TTr]. Let OpG(7rr) denote the sum of the objectives of the optimal solutions to the single-sink 
and single-source problems, and let Opt(n) := Ylr&v Optr(vrr). Let Optj^^g denote the minimum Opt(n) 

*The word tree is misleading since the final solution need not be a tree in directed graphs. Nevertheless, we continue using this 
term for historical reasons. Junction trees were originally proposed for undirected graphs, where the solution is indeed a tree. 

®The single-source problem in a directed graph is identical to the single-sink problem with all the edges reversed in direction. 
For undirected graphs, both problems are on the same graph. 
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over all partitions. We eall this solution the optimum junction-tree solution for this instanee.^*^ Clearly, 
Optjunc > Opt. The junction-tree approximation factor of X is defined to be the ratio Optj^^^/Opt. 

LP Relaxation. We now deseribe a natural flow-based LP relaxafion for fhe single-sink buy-af-bulk problem 
for an insfanee X = (G, T) where X is a sef of terminals fhaf need fo be eonneefed fo fhe roof r. 


minimize E CeXe + EE«w (SS-BaB LP) 

eSS i e£E 

s.f {fi{e): e G E{G)} defines a flow from Si fo r of value 1 Vsj G T 
fi{e) <Xe Ve G X 

Xe > 0, fi{e) >0 ye e E 

Reeall fhaf fhe integrality gap of (SS-BaB LP) on fhe insfanee X = (G, T) is defined fo be fhe rafio of Opt 
to the optimal value of the LP (SS-BaB LP). Also, we define the integrality gap for the graph G to be the 
worst ease integrality gap (over all requests T on graph G) of the eorresponding insfanee X = (G, T). 

2.3 Our Results 

Main Technical Theorem and its Applications. Now we are ready to state our main theorem; the proof is 
in Seetion 3. We say that an online algorithm is y-eompetitive for a graph G if, for any sequenee of requests 
X, the online algorithm for buy-at-bulk returns a solution within a y-faetor of Opt(X), where X = (G, A). 

Theorem 2 (Reduction to Single-Sink Online Algorithms). Fix an instance X = (G,X) of the MC-BB 
problem. Suppose the following three conditions hold. 

(i) The junction-tree approximation factor of X is at most a. 

(ii) The integrality gap of (SS-BaB LP) on any single-sink instance on graph G is at most /3. 

( in) There is a j-competitive online SS-BB algorithm for any instance on graph G that runs in time T. 
Then there is an online algorithm for X running in time poly{n,T) whose competitive ratio is 0{al3'y ■ 
polylog(n)). 

Using this theorem, we ean immediately obtain the following new results mentioned in the introduetion. 

Theorem 3 (Undireeted Edge-weighted Buy-at-Bulk), There is a po\y\og{n)-competitive, polynomial time 
randomized online algorithm for the MC-E-BB problem. 

Proof: The theorem follows direetly by eombining Theorem 2 with the following results from previous 
work. Chekuri et al. [15] prove that the junetion-tree approximation faetor for the undireeted edge-weighted 
buy-at-bulk problem is 0(log k). Chekuri et al. [16] prove that the integrality gap of (SS-BaB LP) in undi¬ 
reeted edge-weighted graphs is 0(log k). Meyerson [28] gives a randomized polynomial time online algo¬ 
rithm for the single-sink buy-at-bulk problem with eompetitive ratio 0(log^ n). ■ 

Theorem 4 (Undireeted Node-weighted Buy-at-Bulk). For any constant e > 0, there is an 0{k^pol'ylog{n))- 
competitive, randomized online algorithm for MC-N-BB with running time As a corollary, this 

yields a polylog{n)-competitive, quasi-polynomial time algorithm for this problem. 

'®Note that copies of the same edge appearing in multiple single-sink solutions are treated as distinct edges in the junction-tree 
solution. Hence, decomposing the optimal multicommodity solution into its constituent paths does not yield Optj,,^,. = Opt. 
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Theorem 5 (Directed Buy-at-Bulk), For any constant e > 0, there is an 0[k^^^^'^polY\og{n)))-competitive, 
polynomial time online algorithm for the MC-D-BB. 

We again use Theorem 2 to prove the above theorems. However, unlike for MC-E-BB, we are not aware of 
any online algorithms for the SS-N-BB and SS-D-BB problems. We therefore first give online algorithms 
for these problems, and then use Theorem 2; the details appear in Sections 5 and 4. 

Finally, we can almost directly use Theorem 2 to also obtain matching results for prize-collecting versions 
of the above problems. Recall that in a prize-collecting problem, every terminal-pair also comes with a 
penalty qi, and the algorithm can opt to not satisfy the request by incurring this value in the objective. We 
give the extension of our results to the corresponding prize-collecting problems in Section 6 . 

Theorem 6 . For each of the above problems, there is an online algorithm with matching running time and 
competitive ratio for the corresponding prize-collecting version. 

In addition to the new results mentioned above, we can also use Theorem 2 to give alternative proofs 
(with slightly worse polylog factors) of some recent results in online network design. By combining Theo¬ 
rem 2 with the polylog (n)-competitive algorithm for online group Steiner Tree due to Alon et al. [2], we ob¬ 
tain a polylog(n)-competitive online algorithm for the group Steiner forest problem - a result shown earlier 
by Naor et al. [31]. Similarly, by combining Theorem 2 with the polylog(n)-competitive online algorithm 
for the node-weighted Steiner tree problem due to Naor et al. [31], we obtain a polylog(n)-competitive 
online algorithm for the node-weighted Steiner forest problem - a result shown earlier by Hajiaghayi et al. 
[24]. 

Height Reduction Theorem. One of the technical tools that we use repeatedly in this paper is the following 
result, which builds on the work of Helvig et al. [25]. We give the proof in Appendix A. 

Theorem 7. Given a directed graph G = {V, E) with edge costs Ce and lengths for all h > 0, we can 
efficiently find an upward directed, layered graph G'^ on {h + 1 ) levels and edges (with new costs and 
lengths) only between successive levels going from bottom (level h) to top (level 0 ), such that each layer has 
n vertices corresponding to the vertices of G, and, for any set of terminals X and any root vertex r, 

(i) the optimal objective value of the single-sink buy-at-bulk problem to connect X (at level h) with 
r (at level 0) on the graph G'^ is at most 0(hk}/^)(j), where f is the objective value of an optimal 
solution of the same instance on the original graph G; 

(ii) given a integral (resp. fractional solution) of objective value f for the single-sink buy-at-bulk 
problem to connect X with r on the graph G'^, we can efficiently recover an integral (resp. fractional 
solution) of objective value at most ffor the problem on the original graph G. 

Likewise, we can obtain a downward directed, layered graph on (h l)-levels with edges going from 

top to bottom, with the same properties as above except for single-source instances instead. 

3 Proof of Theorem 2 (Online Reduction to Single-Sink Instances) 

There are three main steps in the proof. In Section 3.1, we describe the composite LP which is a relaxation 
of optimal junction-tree solutions (for technical reasons, we first need to pre-process the graph). Next, in 
Section 3.2, we show how to fractionally solve the LP online. Third, in Section 3.3, we show how to partially 
round the LP online. The resulting solution then decomposes as fractional solutions to different single-sink 
instances. Finally, we use the bounded integrality gap and the online algorithm for SS-BB to wrap up the 
proof in Section 3.4. 



3.1 The Composite-LP Relaxation: MC-BaB LP 

We first apply Theorem 7 with h = 0(log n) to obtain layered graphs (resp., of height 0(log n) 

where all the edges are direeted upward (resp. downward); see Figure 3 for an illustration. The reason for 
this preproeessing is that the length of the (s*, U) paths appear as a faetor in our final eompetitive ratio and 
the above step bounds it to a logarithmie faetor. Reeall that the graph (resp., G'^°®“) approximately 
preserves the single-sink (resp., single-souree) solutions for any set of terminals and any root. After this 
step, we ean imagine that all the roots (of the single-sink instanees we will solve) are vertiees in level 0 , 
and all the terminals will be vertiees in level h = 0(logn). For elarity of presentation, we refer to the root 
and terminal vertiees by the same name in both G'^p and (even though the graphs are eompletely 

disjoint). Overloading notation, let V denote the vertiees in level 0 in both G“p and O'down^ l-j^g 

union of the edge sets of G“p and Furthermore, the eost Cg and length £e of these edges are inherited 

from Theorem 7. 



GPP 


(terminals arrive here) 


(^down 


Figure 3: The graphs G“p and with h = 3, where the original graph G has 6 vertiees {vi,V 2 , ■ ■ ■ iVq}. 

Now, using the junetion-tree deeomposition with approximation faetor a, we get the following lemma. 

Lemma 8. There exists a set R* C V of root vertices, a partition If* := {vr^ : r £ R*} of the terminal- 
pairs in X, a collection of in-trees : r £ R*} rooted at r in G“p, and a collection of out-trees 

{jidown . 7 - g rooted at r in G^°’"“ such that 

(i) Each {si,ti) £ X belongs to iTrfor some r £ R*. 

(ii) For each r £ R*, the in-tree is a feasible solution to the single-sink buy-at-bulk problem con¬ 
necting {sj : {si,ti) £ TTr} to r in G'^P; likewise, the out-tree A a feasible solution to the 

single-source buy-at-bulk problem connecting r to {f : (s*, ti) £ tt^} in 

(Hi) The sum of objective values of the single-source and single-sink solutions is at most O(alogn) • 
Opt(X). 

Proof: Sinee the junetion-tree approximation ratio of the given instanee is a, there exists a junetion-tree 
solution given by a set of roots R* and a partition Ft* := {vr^ : r £ R*] sueh that the total objeetive value 
of all the single-sink and single-souree junetion-trees is at most aOpt(X). Moreover, by Theorem 7, beeause 
we ehoose the height h = 0(logn), the objeetive value of eaeh single-sink and single-souree solution in G 
inereases by a faetor of O(logn) in G'^p and G'^°"“ respeetively. Thus, the overall objeetive value of the 
resulting junetion-trees in the graphs G“p and is 0{a log n) ■ Opt(X). ■ 
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The above lemma motivates the LP relaxation given in Fig. 4 whieh seeks to assign each {si,ti) pair to 
some rooted instance, and then minimizes the total fractional objective value of the rooted instances. Each 
individual rooted instance is represented by an inner-LP (see the boxed constraints in Fig. 4). 

minimize ^ ^ Cexl + ^ + f^u)) (MC-BaB LP) 

r&V e&E reReeE 

s.t ^ Zir >1 Vi (2a) 

rev 

Zir > 0 (2b) 

(2c) 
( 2 d) 
( 2 e) 
( 2 f) 
( 2 g) 


Figure 4: Composite LP for MC-BB. Equations (2a) and (2b) form the outer-LP; (2c)-(2g) form the inner- 
EPs. 

In the LP, Zir denotes the extent to which the pair {si, U) chooses root r to route its flow. Within each 
inner-LP corresponding to a root r, {x^} are the variables which denotes whether edge e is used to route flow 
in the corresponding rooted instance, and fj , (resp. fT , 0 denotes the amount of flow Si sends (resp., 
ti receives) along e to (resp., from) root r. Observe that if the Zir variables are integral, then the inner-LP 
corresponding to every root r constitutes a feasible fractional solution to (SS-BaB LP) for the single-sink 
instance X' = (G“P, X') where X' = {(sj, r): Zir = 1} and the single-source instance X" = (O^own^ 
where X" = {{r,ti) ■ Zir = 1}- The next lemma, which bounds the optimal value of the MC-BaB LP, 
follows directly from Lemma 8 . 

Lemma 9. The optimum value of (MC-BaB LP) is 0{a\ogn) • Opt(X). 

3.2 An online fractional algorithm for the MC-BaB LP 

Theorem 10. There is a randomized, polynomial-time online algorithm that returns a feasible fractional 
solution for (MC-BaB LP) of value at most 0{a log^ n) • Opt(X). 

In the remainder of the subsection, we prove the above theorem. We remark that the overall reduction 
uses Theorem 10 as a black-box and the time-constrained reader can skip the proof and move to Section 3.3. 

To simplify the exposition, we assume that we know the cost of an optimal solution Opt(X) up to a 
constant factor, using a standard doubling trick'^ Once we know Opt(X), by re-scaling all the parameters 

"Suppose our online algorithm has a competitive ratio of a, and the true cost of an optimal solution is c*. Then, we begin 
with an initial guess for the optimal cost, and run the online algorithm assuming this guess is the correct estimate for c*. If our 
online algorithm fails to find a feasible solution of cost at most a times the current guess, we double our guess and run the online 
algorithm again. Eventually, our guess will exceed the optimal cost c* by at most a factor of two, and for this guess, the algorithm 


{/(e define a flow from s* to r of value Zir in 

Mi, Mr G V 

{/(e ti)} define a flow from r to U of value Zir in (7''°'^“ 

Mi, Mr G V 

fks.) < < 

Mi,Me,Mr G V 

fUu) < < 

Mi,Me,Mr G V 

< > 0, /f) > 0 
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in the problem, we may assume that it equals 1. Next, we delete any edge in or G'iown 
Ce or length larger than 1 as sueh edges eannot partieipate in any optimal solution. Subsequently, we 
initialize all Xg variables to 1/n^. Likewise, we initialize all Zir variables to 1/n® and also send an initial 
flow of 1/n® from eaeh Si to r in on an arbitrary flow path from Si to r and likewise from r to U in 
(^down setting ensures that the eost of the initial solution is o(l). 

In the following, we partition the edge set E into disjoint sets {Ej : 0 < j < h — 1}, where Ej denotes 
the set of edges in E between levels j and j + 1. Furthermore, for elarity of exposition, we deseribe a 
‘eontinuous-time’ version of the algorithm where we inerease the variables as a funetion of time. We note 
that this algorithm ean easily be diseretized for a polynomial-time^^ implementation. The algorithm is given 
as Algorithm 1. 


Algorithm 1 Online Fraetional Algorithm for (MC-BaB LP) 

When a terminal-pair (s*, L) arrives, we update the LP solution using the following steps: 

(1) Let Ri denote the set of roots r in level 0 sueh that s* is eonneeted to r in and r is eonneeted to ti 
in Gdown For eaeh r G Ri, initialize a flow of value 1 /n® using any arbitrary flow path from Si to r in 

and likewise from r to U in Also set Zir = 1/n^ for these roots. 

(2) Repeat the following while ^ir < I- 

(a) Call an edge e G E tight for root r if Xg = /(^g ^ or Xg = ^ 

(b) Edge Update: For all tight edges e e E, update Xg at the rate := 

(e) Flow Update: Solve the following min-eost max-flow problem for eaeh r ^ Ri\ maximize A sueh 
that 

- there exists a flow {^l^g sending A units of flow from Si to r in 

- there exists a flow {^l^g sending A units of flow from r to ti in G‘^°'^“^, 

- Capacity constraints'. < TheAlgorithm(Algorithm l).x^/ce and < x^Jce 

for all tight edges e, 

- Cost constraint: Yje ' 9\e,Si) - Zle 

(d) Update at the rate ■= 9\e^siy fie.u) update 

Zir at the rate = A. 


We inerease the x variables on tight edges at a rate inversely proportional to their eost, similar to the well- 
known online set eover algorithm [3]. However, the “flow eonstraints” are not pure paeking (or eovering) 
eonstraints and there is no general-purpose way of handling them. Indeed, we determine the rate of inerease 
of the flow variables by solving an auxiliary min-eost max-flow subroutine whieh routes ineremental flows of 
equal value from Si to r in G'^p and from r to ti in respeeting eapaeity eonstraints (i.e., for edges that 

are tight, the ineremental flow is at most the rate of inerease of x). This maintains feasibility in the inner LP. 
Moreover, to bound the rate of inerease in objeetive, we enforee that the total length of the ineremental flow 

will compute a feasible solution of cost at most 2ac*. Moreover, since our guesses double every time, the total cost of the edges 
bought by the online algorithm over all the runs across different guesses is at most 0(ac*). 

'^The polynomial is in the size of the input to this algorithm, which for some of our algorithms/results is quasi-polynomial in the 
size of the actual problem instance as stated in the introduction. 
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is at most Zir (this is the “cost” constraint in the min-cost max-flow problem). We stress that the incremental 
flows from the auxiliary problem dictate the rate at which we increase the original flow variables in the LP. 
The final solution is feasible since the algorithm runs until the outer-LP constraint is satisfied. 

First, note that the total cost of initialization is o(l) over all the edge and flow variables. So it suffices 
to bound the cost of the updates. The next lemma relates the total cost of the updates to the total time r for 
which the algorithm runs, and the subsequent lemma bounds r in terms of Opt(X). 

Lemma 11. The LP objective value at the end of the above algorithm is O(logn) • r, where r is the 
(continuous) time for which the algorithm runs. 

Proof: We show that at any time t, the rate of the increase of the LP objective value in the algorithm is at 
most 0(log n); this proves the lemma.The objective increases because of increase in x variables and flow 
variables /; we bound these separately. 

We first upper bound the objective increase due to the changes in the x variables. Fix a level j and let 
denote the set of tight edges in Ej at time t. By definition, (resp., x^) 

equals the total flow on these edges for the pair (si, ti). Since the edges in Ej form a cut separating s* from 
ti, the total flow across this cut is at most Zir- Since Zir < 1, we have ^ •^he 

rate of increase of each such tight edge is precisely X^/Ce, which implies that the total rate of increase of the 
LP value due to the increase of x is at most 


E E E 

Summing over all levels gives the desired 0(log n) bound. 

Next, we upper bound the objective increase due to the changes in the / variables. When these variables 
are updated, the total rate of increase of the objective due to the lengths of the (si, r) and (r, U) flow paths 
is at most Zir — this is precisely the “cost” constraint in the auxiliary flow problem. Hence the total rate of 
increase of flow lengths is at most 2, completing the proof. ■ 

Given the above lemma, we are left to relate r to Opt(X) in order to complete the proof of Theorem 10. 

Lemma 12. The time duration r of the above algorithm satisfies r = 0{a log^ n) • Opt(X). 

We will need several new definitions and auxiliary lemmas in order to prove Lemma 12. Recall from Lemma 8 
that we can assume that the solution that we are comparing against is the set of junction-trees defined by 
and for r G R*. Also, recall that the terminal-pairs are partitioned by the groups n* = {tt,. :r^R*}. 

For every (si, f) G X, let Pf denote the path from Si to the root r in such that (s*, f) G vr^. Similarly, 
let Pj* denote the path from r to f in Let i{P) = Yhe&p Lemma 8 asserts that 

El E E ) =O(alogn)-0 pt(X) (3) 

r&R* \eer“PuTdown ) 

To bound r against the optimal junction-tree solution, we use two sets of charging clocks: 

'^We remark that the “log n” corresponds to the number of levels in Gh justifying the preprocessing step before the LP descrip¬ 
tion. 
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• We maintain an edge clock on every (e, r) pair sueh that e G r"P or e G i.e., if e is used by 

the optimal junetion-tree solution in the single-souree (or single-sink) instanee eorresponding to r. In 
partieular, note that if an edge e is in multiple junetion-trees, then it has a separate eloek for eaeh sueh 
tree. 

• We maintain a terminal clock on every terminal-pair {si, U) G Af. 

The erueial invariant that we maintain is the following: at any time instant t, at least one clock “ticks,” i.e., 
augments its counter at unit rate. The overall goal would then be to bound the total time for whieh all the 
eharging eloeks ean eumulatively tiek. 

First, we deseribe the rule for the tieking of the eloeks. Fix a time t, and let the terminal-pair {si,ti) be 
the pair that is aetive at time t. Let r denote the root vertex whieh {si,ti) has been assigned to in the optimal 
junetion-tree solution from Lemma 8, i.e., {si, U) G itr. Now, eonsider the flow-paths P*. in and P^_ in 

(^down j^^^g gj^g gj: |.,^^g situations: 

- If any variable Xg is tight for any edge e G P* U P^ at time t, then the edge eloek on the pair (e, r) 
tieks at time t. If there are multiple sueh edges, then all the eorresponding eloeks tiek. 

- Otherwise, both paths are free of tight edges. In this ease, the terminal eloek for (s*, ti) tieks at time 
t. 


Lemma 13. For any pair (e, r) such that e G U its edge clock ticks for 0(ce log n) time. 

Proof: Notiee that Xg is initialized to l/of for all roots r, and inereases at the rate 


dt Cg 


(4) 


at all times when the edge clock on (e, r) ticks. To see why, eonsider a time t when the eloek on (e, r) tieks, 
and let (sj, fj) denote the aetive terminal-pair at time t. It must be that (i) (s*, U) has been assigned to root 
r in n*, and (ii) either Xg = ^ ^ or Xg = ^ y But in this ease, we inerease sueh variables at rate x^^jce 

in our algorithm (Step (2a)). Therefore, we ean infer that the value of Xg would be 1 after the edge eloek 
on e has tieked for time 0(cg logn). But elearly, e eannot be a tight edge for any subsequent terminal-pair 
(sipi) onee Xg reaehes 1; therefore, the edge eloek on (e, r) tieks for 0(ce logn) time overall. ■ 


Lemma 14. For every terminal-pair {si,ti) connected by the optimal junction-tree solution through the root 
vertex r, the total time for which the terminal clock ticks is at most O(logn) • max(i{Pf),i(P^.)). 

Proof: Reeall that if the terminal eloek for (sj, ti) is tieking at time t, then it must mean that no edge is 
tight on either path Pf or P^. In this ease, we show that the variable Zir inereases at a fast enough rate, 
where r is the root {si,ti) is assigned to in the optimal junetion-tree, i.e., {si,ti) G vr^. We show this by 
exhibiting a feasible solution to the auxiliary LP eonsidered in Step (2b) of the algorithm for root r. Indeed, 
send the flow from Si to r along Pf, and likewise from r to U along P^. Also set the value of A to be 
Zirj max(£(P*,), £(P^*J). Clearly, on the edges of these flow paths, we do not have any eapaeity eonstraints 
sinee no edge is tight. So, the only eonstraints are the eost eonstraints whieh are satisfied by the ehoiee of 
A. Henee, the rate of inerease of Zir is at least 


dZij. ^ 

dt ~ 


Zir 

max{e{Pf),i{P*^) 


(5) 


13 




at all times when the terminal clock on (sj, ti) ticks. This proves the claim, for otherwise the variable Zir 
would have reached 1, and the algorithm would have completed processing (si, fi). ■ 

Since at least one clock ticks at all times, the total time clocked is at least r, the duration of the algorithm. 
Lemma 13 and Lemma 14 imply that 

r<0(logn)j;[ {£{Pl)+i{P*)) 

r&R* \eeT“PuTd°wn 

which together with (3) completes the proof of Lemma 12. Theorem 10 follows from Lemma 11 and 
Lemma 12. 

3.3 Partial Online LP Rounding 

We partially round the fractional solution returned hy Theorem 10 to obtain integral values for only the 
outer-LP variables Zir, i.e., each {si,ti) pair is integrally assigned to a root. The inner-LP variables x and 
/ continue to be fractional but represent unit fractional flow from s* to r and r to ti for the {si, U) pairs 
assigned to r. The partial rounding algorithm is given as Algorithm 2. 


Algorithm 2 Online Partial Rounding Algorithm 

(1) Initialization: Each root chooses a threshold G [1/ 2n, 1/(3 log n)] uniformly at random. 

(2) Partial Rounding: At each time, maintain the scaled solution xl = min(l,x(;/r^), = 

min f 1, . Also set Zir = 1 if Zir > r^. 


Theorem 15. The scaled solution (x, f) component-wise dominates a feasible solution to the outer-LP, and 
the expected objective value of the scaled solution {x, f) is at most 0{alog^ n) ■ Opt(X). Moreover, for 
each (sj, tf), there exists at least one root r such that Zir > 1 with probability at least 1 — l/rv'. 

Proof: Since each root r chooses its threshold Tr independently and uniformly at random from [l/2n, 1/logn], 
the probability that Zir = 1 is at least Zir log n (since Zir = 1 if and only if < Zir). Since this is indepen¬ 
dent for different roots, a standard Chernoff-Hoeffding bound application (see, e.g., [30]) shows that each 
{si,ti) pair has Zir = 1 for some root r with probability at least 1 — 1/n^. Moreover, the expected value of 
any variable Xg is given by 


Aogn r 

/ — log ndTr < 0(log^ n)Xg. 

J Tr = l/2n p 

A similar argument shows that the expected values of scaled flow variables are also bounded by 0(log^ n) 
times their values in the fractional solution. This shows that the expected objective value of the (x, /) solu¬ 
tion is at most 0(log^ n) times the value of (x, /); by Theorem 10, the latter is at most 0{a log^ n)Opt(X). 
Combining these facts gives us the desired bound on the value of the scaled solution. 

It remains to show that the scaled solution dominates a feasible solution to the LP. To this end, fix 
some root r and let A/ denote the set of {si,ti) pairs for which Zir = 1. We need to show that installing 
capacities of {/(^g ^y} on the edges can support unit flow from Si to r in for all (s*, U) G AV. Suppose 
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for contradiction that there is is a cut Q separating Si from r of capacity strictly smaller than 1. This 
implies that every edge e £ Q must have ^ ^ < r^; otherwise, we would have an edge e with ^ ^ = 
1, which contradicts our assumption on the cut capacity. But then the value of the min-cut is precisely 
(XlesQ J^les )) which must be at least 1 because of the following two observations: (i) we know that 
{fj^es is a feasible flow from s* to r of value Zir and hence it must be that Yle&Q f{es) — 
since Zir = 1, it must be that Zir > Tr- This contradicts the assumption that the cut capacity is strictly 
smaller than 1. A similar argument shows that the variables {/(^g ^,)} can support unit flow from r to U for 
every {si,ti) with Zir = 1. ■ 


3.4 Wrapping up: Invoking the Single-Sink Online Algorithm 

We are now ready to put all the pieces together and present our overall online multicommodity buy-at-bulk 
algorithm as Algorithm 3. SingleSinkAIg is the online algorithm for SS-BB alluded to in point (iii) of the 
statement of Theorem 2. 

Algorithm 3 Online Multicommodity Buy-at-Bulk Algorithm 
when {si,ti) arrives 

(1) update the fractional solution of the composite LP using the algorithm (Algorithm 1, Section 3.2). 

(2) partially round the solution using algorithm in (Fig. 2, Section 3.3). 

(3) if(3r : Zir > 1): send both s* and ti to the instance of SingleSinkAIg with root r. 

(4) else: buy a trivial shortest path between Si and U on the metric (c + £) and route along this path 


Clearly the algorithm produces a feasible solution; so we now argue about the expected objective value. 
Fix an (s*, U) pair. Since the probability that a terminal-pair is not assigned to a root is < 1/n^ (by The¬ 
orem 15), the expected total contribution of such unassigned terminal-pairs is < 1 = Opt(X). For a root 
r, let TTr be the terminal-pairs assigned to r. We know that {x, f) restricted to tt^ dominates a feasible 
solution in (SS-BaB LP). Letting LPr denote the contribution of this restriction to the overall LP value, 
we get LPr = 0(alog^ n) • Opt(X). By the integrality gap condition, we get that Opt^, i.e. the in¬ 
tegral optimum objective value of the instance generated by r and vr,., is at most /3 • LPr. (Here we are 
using the fact from Theorem 7 that moving to the layered instance does not increase the integrality gap.) 
The objective value of the solution produced by SingleSinkAIg is at most 7 • Opt^, where 7 is the com¬ 
petitive ratio of SingleSinkAIg. Putting these observations together, we conclude that the overall objective 
value of the solution returned by the online algorithm is 0{aj3'^ log^ n) • Opt(X). This completes the proof 
of Theorem 2. 

4 Online Directed Buy-at-Bulk 

In this section, we prove Theorem 5. A natural approach is to use the reduction given by Theorem 2. To this 
end, we need to establish the following: the existence of a junction-tree scheme with a good approximation; 
a good upper bound on the integrality gap for single-sink instances of the LP given in Section 2; and an 
online algorithm for single-sink instances with a good competitive ratio. 
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Extending the work of Chekuri et al. [13], Antonakopoulos [5] shows the existenee of a junetion-tree 
seheme with approximation 0{'/k). Unfortunately, the integrality gap of the LP relaxation is not very well 
understood even for Steiner tree instanees; [35] gives an U(\/fe) lower bound^^ on the integrality gap for the 
Steiner tree problem and no suitable upper bound is known. We overeome this diffieulty as follows. Instead 
of working with general graphs, we pre-proeess the instanee and obtain a tree-like graph for whieh we ean 
show that the LP has a good integrality gap. Finally, we give the first non-trivial online algorithm for the 
directed single-sink buy-at-bulk problem. These results, together with our reduction (Theorem 2), imply the 
online algorithm for MC-D-BB. 

We devote the rest of this section to the proof of Theorem 5; to aid the reader, we restate the theorem 
below. 

Theorem 16. For any constant e > 0, there is a 0{k^~^^]yo\'ylog{n))-competitive, polynomial time ran¬ 
domized online algorithm for the general buy-at-bulk problem. 

Pre-processing step. We first give our reduction from general instances of the directed buy-at-bulk problem 
to much more structured instances; the reduction loses a factor of 0(/c 2 +^) in the approximation ratio. 

Let h = [1/e]. Given an instance X = {G,X) of the directed buy-at-bulk problem, we map it to a 
tree-like instance J' = {H,X) as follows. We start by applying Theorem 7 to G to obtain the graphs 
and recall that these graphs are layered {h -\- 1)-level graphs with n vertices (corresponding to the 

vertices in G) in each level, and the levels are numbered 0,1,..., /i with 0 being called the root level. The 
graph has edges directed from higher numbered levels to lower numbered levels, and has edges 

in the opposite direction. To facilitate the construction of the graph H, we now create n trees from and 
n trees from G^°’"“ as follows. 

For every “root vertex” r at level 0 in G“p (resp. G'^°"“), the tree (resp. is constructed as 

follows: 

• The 0* layer of has just one vertex - the root r. 

• For each i such that 1 < i < h, the f-th layer of contains all (i -|- 1)-length tuples (r, ui,..., n*) 
where vj is a vertex present in the y-th layer of G'^p. 

• For every edge e = {vi, Vi-i) G G'^p, there is an arc from (r, ui,..., Vi-i,Vi) to (r, ui,..., Uj-i) 
inheriting the same cost Ce and length 

Therefore each tree r“P is an in-arborescence, with all edges directed towards the root. The tree is 

constructed analogously except all edges are directed away from the root. In the following, we use the term 
leaves to refer to the vertices on layer h of these trees. 

After performing the above operation for every root vertex r in level 0 of Gf^ and we have 2n 

trees. Then the final graph H is obtained as follows (see Figure 5). For each root r G U, we first add an arc 
from the root of r“P to of zero cost and length. Finally, for every {si, U) pair, we add the vertices Si 

and ti to H and the following arcs connecting them to the trees: for each tree r“P, we add an arc from Si 
to each leaf of of the form {r,vi,, Vh) with Vh = sp, for each tree we add an arc to U from 

each leaf of of the form (r, ui,..., Vh) with Vh = f. These new arcs have zero cost and length (i.e., 
Ce = ie = 0). 

This completes the construction of H. Note that the graph H has vertices and a similar number 
of edges. Our new instance is = {H,X) and we will apply Theorem 2 to this instance. 

We first relate the objective values of J and X. 

'"'However, in these instances, n is exponentially large in k. So, they do not rule out a polylog(n) upper bound. 
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Figure 5: Construction of graph H. 


Lemma 17. Every feasible solution for J is a junction-tree solution. 

Proof: Note that any {si, ti) path in H has the following structure: Si connects to a leaf node of for 
some r £ V, then continues to the root r, then traverse the edge to the root of then goes down to 

a leaf of and finally connects to f. Thus, for any feasible solution for J', the {si,ti) pairs can be 

partitioned based on the root r through which they connect. ■ 

Lemma 18. Any feasible solution for J can be mapped to a feasible solution — in fact, a junction-tree 
solution —forX of equal or smaller objective value. 

Proof: Note that from the previous lemma, any feasible solution S' in .T” is a junction-tree solution. There¬ 
fore, there is a partition of the (s*, f) pairs depending on which root vertex they are using to connect. 
Moreover, it follows from our construction of the trees in H that any edge in (resp. corresponds 

to an edge in (resp. Therefore, if we map each edge appearing in solution S to its correspond¬ 
ing edge in or we obtain a mapping from each junction tree of S rooted at r to a junction tree in 

G“P U rooted at r that is connecting the same subset of pairs. Finally, by Theorem 7, each junction 

tree in U rooted at r can be mapped, without increasing the objective value, to a junction tree in 

G rooted at r that is connecting the same subset of pairs. This completes the proof of the lemma. ■ 

Lemma 19. Opt(,y) < where is the objective value of an optimal junction- 

tree solution for X. 

Proof: Consider the optimal junction tree solution for X. Let the optimum partition be Ft = (vr^i,.. •, 
where R* = {ri, r 2 ,..., Vq] is the set of roots of the junction trees. For each r G R*, let Xr = {s* : 
{siRi) £ TTr} be the set of sources of and be the corresponding sinks. From Theorem 7, we know 
that the optimum objective value of any one single-sink problem connecting to r in G“p is at most 
0{hk^/^) times the objective value of the optimum solution connecting each source in X^ to r. An analo¬ 
gous upper bound holds for every optimal single-source solution connecting r to each sink in Ij.. Therefore, 
we get that the total sum of objective values of each of the junction trees in G'^p and is at most 
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O(/i.A:^/^)0ptjync(2^)- Now notice that any solution Sq for (C^p, Xr) can easily be “simulated” by a solu¬ 
tion St in the tree indeed, for every root-vertex path (r, ni, ^ 2 ,..., fj) in the solution Sq, include the 
edge from (r, ni, ^ 2 ,..., Wj) to (r, ^ 1 ,^ 2 ,..., Vi-i) in St (recall the vertices in exactly correspond to 
such root-vertex paths). It is easy to see that the objective value of the solution St in is the same as that 
of Sg- Similarly, any solution for (G'down^ simulated in with the same objective value. It 

follows that there is a feasible solution in J of objective value at most 0(/i/c^/^)0ptju„^(X). ■ 

Corollary 20. Opt(J') < O(A:^+^)0pt(X). 

Proof Sketch: Antonakopoulos [5] shows that there exists a junction-tree solution of cost at most O (Vk) 0 pt (X). 

The corollary follows from this work and the fact that we set h = 0(l/e). ■ 

Now we are ready to show that the new instance J' has the properties required by the reduction, i.e.. 
Theorem 2 can be applied. In the following lemma, a single-source (resp. single-sink) sub-instance 
J' = {H,T,v) of J' = {H,X) is a single-source (resp. single-sink) instance of the following form: 
the graph is the same as in J, namely H ; the set of terminals T is a subset of the sources (resp. sinks) of 
the terminals T need to be connected to a root vertex v ^ V {H) \ {s*, U: i ^ [A:]}. 

Lemma 21, Let J be the instance described above. Let a, /?, 7 be as in the statement of Theorem 2. We 
have 

(i) The junction-tree approximation factor of J is 1 (i.e., a = 1). 

(ii) The integrality gap of (SS-BaB LP)/or any single-sink/source sub-instance J' = (H, T, v) is 0(h log n log k) 
(i.e., (3 = 0{h\ogn\ogk)). 

(Hi) There is a 0{h?\o^ n\ogk)-competitive algorithm for any single-sink/source sub-instance J' = 

(H, T, v) (i.e., 7 = 0{h? log^ n log k)). 

Proof: Property (i) follows from Lemma 17. Thus we focus on proving (ii) and (iii). In the following, we 
assume that we are working with a single-sink sub-instance J’' of J\ the proof is very similar for single¬ 
source sub-instances and we omit it. 

In order to show (ii) and (iii), we will map the sub-instance J' = (77, X, v) to an instance of the group 
Steiner tree problem on a tree as follows. Since v £ V (H) \ {sj, f: i £ [fc]}, we have v £ U 
for some r. We first consider the case when v £ In order to define fhe free of fhe group Sfeiner free 
insfance, we sfarf wifh fhe subfree X of roofed af v. We add fhe following ‘dangling’ edges fo Ty for 
each source Si £ T: for each leaf vertex u £ such fhaf Si has an edge in fo u, we add a new vertex 
u' fo Ty and connecf if fo u. Lef T be fhe resulting free. We assign weighfs fo fhe edges of T as follows. Each 
of fhe old edges e £ Ty receives a weighl equal fo Cg. Each of fhe new edges uu' £ E(T)\ E(Ty) receives a 
weighl equal fo i{Pu), where Py is fhe pafh of from u fo v. Einally, we define fhe following groups: for 
each source Sj G T, we infroduce a group Si consisfing of all fhe new verfices u' £ V(T)\V (Ty ) such fhaf 
Si is connecfed fo ifs parfner u in (fhaf is, has an edge from Si fo u). In fhe resulting group Sfeiner 
free insfance, fhe goal is fo connecf all of fhe groups 5 = {5^: Sj G T} fo fhe roof v using a minimum 
weighf subfree of T; we lef (T, S, v) denote fhis insfance. 

Now fhe key claim is fhaf fhe feasible solutions fo fhe single-sink buy-af-bulk (77, T, v) are in a one- 
fo-one correspondence wifh feasible solution fo fhe group Sfeiner free insfance {T,S,v)', moreover, fhe 
objecfive value of a solufion fo fhe former is equal fo fhe weighf of fhe solution fo fhe laffer. To see fhis, 
consider a feasible solufion V for fhe single-sink buy-af-bulk insfance. Note fhaf, for each source Si £ T,V 
has a pafh connecfing s* fo n; if follows from our consfrucfion of 77 fhaf fhis pafh consisfs of an edge from 
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Si to a leaf u of followed by the unique path in from u to v. Thus we ean eonstruet a feasible group 
Steiner tree solution by eonneeting eaeh group Si using the path of T from u' to v, where u' is the partner of 
the leaf u of through whieh Si eonneets to the root in the buy-at-bulk solution. The weight of the edge 
u'u eaptures the £-eost of Sj’s path and the weight of the path from uto v eaptures the c-eost of Sj’s path. 

Moreover, we ean apply the same argument to fraetional solutions to the two problems and show that 
there is a bijeetion between feasible fraetional solutions to (SS-BaB LP) and feasible fraetional solutions 
to the LP relaxation for group Steiner tree of Garg, Konjevod, and Ravi [19]; as before, these eorre- 
sponding solutions have the same objeetive values. Now the desired upper bound on the integrality gap 
of (SS-BaB LP) follows from the work of [19] who showed that the integrality gap of the group Steiner 
tree LP is 0{logNlogK) where N = maxj \Si\ is the maximum size of a group and K is the number of 
groups. In our setting, K < |Tf| < k and N < n^. Therefore the integrality gap is 0{h\ogn\ogk), whieh 
establishes property (ii). 

Moreover, notiee that the above reduetion ean also be used to obtain an online algorithm for the single¬ 
sink (and single-souree) sub-instanees. Indeed, we simply use the online group Steiner tree algorithm of 
Alon et al. [1] whieh has a eompetitive ratio 0(log^ NlogK) = 0{h? log^ nlog k). This proves property 
(iii). ■ 

Now we are ready to eomplete the proof of Theorem 16. 

Proof of Theorem 16: Given the instanee X = (G, Tf), we eonstruet the graph H as deseribed above; the 
time taken to do so is We pass the instanee {H, X) to Theorem 2 and, using Lemma 21, we obtain 

an online algorithm that, for any eolleetion of pairs X and any adversarial ordering of X, returns a solution 
of eost 0(polylog(n)) • Opt(J"). By Lemma 19 and Corollary 20, we ean map solutions for {H,X) to 
solutions for (G, Tf). ■ 

We note that the approaeh deseribed above also gives us new online algorithms for the single-sink buy-at- 
bulk problem on direeted graphs. For single-sink instanees, the junetion-tree approximation is equal to 1 
and thus we save a faetor of s/k. 

Corollary 22. There is a polynomial time 0{k^ ■ polylog{n))-competitive online algorithm for the single¬ 
sink (or single-source) buy-at-bulk problem on directed graphs. The competitive ratio can be improved to 
polylog (n) if the running time can be quasi-polynomial in n. 


5 Online Single-sink, Undirected, Node-weighted Buy-at-Bulk 

In this seetion, we prove Theorem 4. Again, we use our main theorem (Theorem 2) and reduee the multi- 
eommodity buy-at-bulk problem to the single-sink version of the problem. We eombine the reduetion theo¬ 
rem with the following results from previous work. Chekuri et al. [13] show the existenee of a junetion-tree 
seheme with approximation faetor 0{logk). Moreover, the natural LP relaxation for the single-sink buy- 
at-bulk problem on graphs with node eosts is also 0(log/c) [15]. For the single-souree online algorithm, 
we resort to the algorithms from the previous seetion for the more general direeted single-sink buy-at-bulk 
problem (see Corollary 22). We now obtain the desired result by the following parameter settings: 

a = 0{logk) 

(3 = Oilogk) 

7 = 0(polylog(n)) 

T = 
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6 Online Prize-Collecting Buy-at-Bulk 


In the prize-collecting version of the buy-at-bulk problem, each terminal pair {si,ti) also comes with a 
penalty Qi and the algorithm may choose not to serve this request and incur the penalty in the total cost. We 
show that our online reduction framework (Theorem 2) can be easily modified to handle prize-collecting 
versions as follows. 

Theorem 23. Let X be a buy-at-bulk instance and suppose the three conditions of Theorem 2 hold. Then 
there is an 0{aj3y • y)o\y\og{n))-competitive online algorithm for the online, prize-collecting buy-at-bulk 
problem on X with arbitrary penalties. 

Proof: We closely follow the proof of Theorem 2. The first difference is in the LP-formulation. Now, for 
each (sj, f) pair we have an extra variable Zi^ which indicates whether we choose to discard this pair (and 
pay the corresponding penalty) or not. We point out the differences with (MC-BaB LP). The new objective 
function is 

minimize ^ +/(m,)) + 

r&V e&E isi,ti)£X r&Re&E 

and (2a) is replaced by 

T ^i,0 

rev 

Observe that the optimum value of this modified LP is af mosf 0 (q; log n) fimes fhe opfimum: sef Zj o = 1 
for fhe pairs fhe infegral opfimum solufion does nof connecf, and for fhe resf apply Lemma 9. Also observe 
fhaf fhe modified LP can be fhoughf of fhe old LP on a modified insfance where fhe graphs and 
(obfained from Theorem 7) have anofher verfex “0” af fhe roof level, and each Sj has a direcf pafh from s* fo 
“0” in wifh fofal lengfh qi/2 and no fixed cosf, and similarly, each f has a pafh from “0” fo U in 
wifh fofal lengfh qi/2 and no fixed cosf. The resf of fhe proof now follows exacfly as in Secfion 3, by also 
including fhe special verfex as a possible roof while rounding fo make fhe oufer LP variables infegral. ■ 

7 Conclusion 

In fhis paper, we gave fhe firsl polylogarilhmic-compelifive online algorifhms for fhe non-uniform mulfi- 
commodify buy-af-bulk problem. Our resulf is a corollary of a generic online reduction technique fhaf we 
proposed in fhis paper for converting a mulficommodify insfance info several single-sink insfances, which 
are often easier fo design algorifhms for. We believe fhaf fhis reduction will have ofher applicafions be¬ 
yond fhe buy-af-bulk framework, and illusfrafe fhis by showing fhaf recenf resulfs on online node-weighfed 
Steiner foresf and online generalized connecfivify direcfly follow from our reducfion fheorem. Our work 
also opens up new directions for fufure research. For insfance, our algorifhm for fhe node-weighfed problem 
runs in quasi-polynomial time, and a concrefe open quesfion is fo gef a polynomial-time polylogarifhmic- 
compefifive algorifhm for fhe SS-N-BB problem (fhis suffices for MC-N-BB as well by our main fheorem). 
Anofher technical question concerns non-uniform demands. While our algorifhm can be extended fo fhe 
case of non-uniform demands, fhe approximafion ratio incurs an addifional O (log 79) factor, where D is 
fhe rafio of fhe largesf fo fhe smallesf demand. If would be inferesfing fo eliminate fhis dependence on 79 
since fhe corresponding offline resulfs do nof have fhis dependence. More generally, a broader quesfion is 
fo investigate ofher mixed packing-covering LPs fhaf can be solved and rounded online. 


20 



Acknowledgements 


D. Panigrahi is supported in part by NSF Award CCF-1527084, a Google Faeulty Researeh Award, and a 
Yahoo FREP Award. 


References 

[1] N. Alon, B. Awerbueh, Y. Azar, N. Buehbinder, and J. Naor. A general approaeh to online network op¬ 
timization problems. In Proceedings, ACM-SIAM Symposium on Discrete Algorithms (SODA), pages 
577-586, 2004. 

[2] N. Alon, B. Awerbueh, Y. Azar, N. Buehbinder, and J. Naor. A general approaeh to online network 
optimization problems. ACM Trans. onAlg., 2(4):640-660, 2006. 

[3] N. Alon, B. Awerbueh, Y. Azar, N. Buehbinder, and J. Naor. The online set eover problem. SIAM J. 
Comput, 39(2):361-370, 2009. 

[4] M. Andrews. Hardness of buy-at-bulk network design. In Proceedings, IEEE Symposium on Eounda- 
tions of Computer Science (EOCS), pages 115-124, 2004. 

[5] S. Antonakopoulos. Approximating direeted buy-at-bulk network design. In Proceedings, Workshop 
on Approximation and Online Algorithms (WAOA), pages 13-24, 2010. 

[6] B. Awerbueh and Y. Azar. Buy-at-bulk network design. In Proceedings, IEEE Symposium on Eounda- 
tions of Computer Science (EOCS), pages 542-547, 1997. 

[7] Y. Azar, U. Bhaskar, L. Fleiseher, and D. Panigrahi. Online mixed paeking and eovering. In Proceed¬ 
ings, ACM-SIAM Symposium on Discrete Algorithms (SODA), pages 85-100, 2013. 

[8] P. Berman, A. Bhattaeharyya, K. Makaryehev, S. Raskhodnikova, and G. Yaroslavtsev. Approximation 
algorithms for spanner problems and direeted Steiner forest. Inform, and Comput., 222:93-107, 2013. 

[9] P. Berman and C. Coulston. On-line algorithms for Steiner tree problems (extended abstraet). In 
Proceedings, ACM Symp. on Theory of Computing (STOC), pages 344—353, 1997. 

[10] N. Buehbinder and J. Naor. Online primal-dual algorithms for eovering and paeking. Math. Open Res., 
34(2):270-286, 2009. 

[11] M. Charikar, C. Chekuri, T. Cheung, Z. Dai, A. Goel, S. Cuba, and M. Li. Approximation algorithms 
for direeted Steiner problems. J. Algorithms, 33(1):73-91, 1999. 

[12] M. Charikar and A. Karagiozova. On non-uniform multieommodity buy-at-bulk network design. In 
Proceedings, ACM Symp. on Theory of Computing (STOC), pages 176-182, 2005. 

[13] C. Chekuri, G. Even, A. Gupta, and D. Segev. Set eonneetivity problems in undireeted graphs and the 
direeted Steiner network problem. ACM Transactions on Algorithms, 7(2): 18, 2011. 

[14] C. Chekuri, M. T. Hajiaghayi, G. Kortsarz, and M. R. Salavatipour. Approximation algorithms for 
node-weighted buy-at-bulk network design. In Proceedings, ACM-SIAM Symposium on Discrete Al¬ 
gorithms (SODA), pages 1265-1274, 2007. 


21 



[15] C. Chekuri, M. T. Hajiaghayi, G. Kortsarz, and M. R. Salavatipour. Approximation algorithms for 
nonuniform buy-at-bulk network design. SIAM J. Comput., 39(5): 1772-1798, 2010. 

[16] C. Chekuri, S. Khanna, and J. Naor. A deterministie algorithm for the eost-distanee problem. In Pro¬ 
ceedings, ACM-SIAM Symposium on Discrete Algorithms (SODA), volume 7, pages 232-233, 2001. 

[17] Y. Dodis and S. Khanna. Design networks with bounded pairwise distanee. Proceedings, ACM Symp. 
on Theory of Computing (STOC), pages 750 - 759, 1999. 

[18] M. Feldman, G. Kortsarz, and Z. Nutov. Improved approximation algorithms for direeted Steiner forest. 
J. Comput. System Sci., 78(l):279-292, 2012. 

[19] N. Garg, G. Konjevod, and R. Ravi. A polylogarithmie approximation algorithm for the group Steiner 
tree problem. J. Algorithms, 37(l):66-84, 2000. 

[20] S. Cuba, A. Meyerson, and K. Munagala. A eonstant faetor approximation for the single sink edge 
installation problem. SIAM J. Comput., 3i{6):2426-2442, 2009. 

[21] A. Gupta, A. Kumar, M. Pal, and T. Roughgarden. Approximation via eost-sharing: A simple approx¬ 
imation algorithm for the multieommodity rent-or-buy problem. In Proceedings, IEEE Symposium on 
Eoundations of Computer Science (EOCS), pages 606-615, 2003. 

[22] A. Gupta, A. Kumar, and T. Roughgarden. Simpler and better approximation algorithms for network 
design. In Proceedings, ACM Symp. on Theory of Computing (STOC), pages 365-372, 2003. 

[23] M. Hajiaghayi, V. Liaghat, and D. Panigrahi. Near-optimal online algorithms for prize-eolleeting 
Steiner problems. In Proceedings, International Colloquium on Automata, Languages and Processing 
(ICALP), pages 576-587, 2014. 

[24] M. T. Hajiaghayi, V. Liaghat, and D. Panigrahi. Online node-weighted Steiner forest and extensions 
via disk paintings. In Proceedings, IEEE Symposium on Eoundations of Computer Science (EOCS), 
pages 558-567, 2013. 

[25] C. S. Helvig, G. Robins, and A. Zelikovsky. An improved approximation seheme for the group Steiner 
problem. Networks, 37(1):8-20, 2001. 

[26] M. Imase and B. M. Waxman. Dynamie Steiner tree problem. SIAM J. Discrete Math., 4(3):369-384, 
1991. 

[27] S. Korman. On the use of randomization in the online set eover problem. M.S. thesis, Weizmann 
Institute of Science, 2005. 

[28] A. Meyerson. Online algorithms for network design. In Proceedings, ACM Symposium on Parallelism 
in Algorithms and Architectures (SPAA), pages 275-280, 2004. 

[29] A. Meyerson, K. Munagala, and S. A. Plotkin. Cost-distanee: Two metrie network design. SIAM J. 
Comput., 38(4): 1648-1659, 2008. 

[30] R. Motwani and P. Raghavan. Randomized Algorithms. Cambridge University Press, 1997. 

[31] J. Naor, D. Panigrahi, and M. Singh. Online node-weighted Steiner tree and related problems. In 
Proceedings, IEEE Symposium on Eoundations of Computer Science (EOCS), pages 210-219, 2011. 


22 



[32] F. S. Salman, J. Cheriyan, R. Ravi, and S. Subramanian. Approximating the single-sink link- 
installation problem in network design. SIAM J. Optimization, 11(3):595-610, 2001. 

[33] K. Talwar. The single-sink buy-at-bulk LP has eonstant integrality gap. In Proceedings, MPS Confer¬ 
ence on Integer Programming and Combinatorial Optimization (IPCO), pages 475-486, 2002. 

[34] A. Zelikovsky. A series of approximation algorithms for the aeyclic direeted Steiner tree problem. 
Algorithmica, 18(1):99-110, 1997. 

[35] L. Zosin and S. Khuller. On directed Steiner trees. In Proceedings, ACM-SIAM Symposium on Discrete 
Algorithms (SODA), pages 59-63, 2002. 

A Reduction to Layered Instances (Proof of Theorem 7) 

In this section, we prove Theorem 7, which is an extension of Zelikovsky’s ‘height reduction lemma’ for the 
buy-at-bulk problem; Zelikovsky’s original Lemma was for a single metric, whereas in our setting there is 
both a cost and a length metric. 

We prove the up-ward case; the down-ward case follows analogously. In order to simplify the notation, 
we remove the superscript For this reduction, we will adapt the notion of layered expansion of a graph, 
which has been in the folklore for many years and has been used recently by several papers (see, e.g.,[13, 
31]). The /i-level layered expansion of G is a layered DAG Gh of h-\-l levels (we index the level 0,1,... ,h) 
defined as follows: 

(i) For each i such that 0 < i < h, the vertices in level i are copies of the vertices of G; we let Vi to 
denote the copy of vertex v at level i. 

(ii) For each i such that 1 < i < h, there is a directed edge from every vertex in level i to every vertex in 

level f — 1. The fixed cosf of an edge (ui, Vi-f) is given by fhaf of fhe shorfesf direcfed pafh from 
Ui fo Vi-i in G according fo fhe mefric Ce -|- The lengfh of fhis edge is sef fo be fhe lengfh 

of fhe pafh in fhe I mefric. 

We now relafe fhe opfimal objecfive values for fhe fwo insfances. One of fhe direcfions of fhe reduction 
is slraighlforward. 

Lemma 24. For any root r and any set of terminals X, if there is a feasible integral/fractional solution 
of objective/LP value f for the single-sink buy-at-bulk problem connecting X to r on the h-level layered 
expansion Gh, then there is a feasible integral/fractional solution of objective/LP at most f for the same 
problem in G. 

Proof: Note fhaf for every edge in Gh, fhere is a corresponding pafh in G wifh fhe properly fhaf fhe sum of 
edge cosls and lenglhs on fhe pafh is al mosl fhe cosf and lengfh of fhe edge in G. Therefore, replacing fhe 
edges in fhe solution for fhe layered graph by fhe corresponding pafhs in G yields a feasible solution in G 
wifhouf increasing fhe overall cosf and lengfh. Nolice fhaf fhe same “embedding” of edges in Gh to paths in 
G can be applied to the fractional solution on Gh as well. This shows property (ii) of the theorem statement. 


The more interesting direction is to show that the optimal objective value on the layered graph Gh can be 
bounded in terms of the optimal objective value on the original graph G. To show this, we will re-purpose 
the so-called “height reduction” lemma of Helvig, Robins, and Zelikovsky [25]. We restate the lemma in a 
form that will be useful for us. 
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Lemma 25. For any in-tree T defined on the edges ofG that is rooted at r and contains all the terminals in 
X, and for any integer h > 1, there is an in-tree T' (on the same vertices as G but over a different edge set) 
that is also rooted at r and contains all the terminals, and has the following properties: 

(i) T' contains h + 1 levels of vertices, i.e., has height h. 

(ii) T' is an k^^^-ary tree, i.e., each non-leaf vertex has children. 

(in) Each edge e' = {u', v') in tree T' corresponds to the unique directed path p^,/ in T from u' to v'. 
Moreover, the number of terminals in the subtree ofT rooted at u' is exactly where e' is an 

edge between levels i — 1 and i ofT. 

(iv) Each edge in T is in at most 2hk^^^ such paths pe' for edges e' G T'. 

For an edge e' G T', suppose we define its eost to be the eost of the path p^', and its length to be the 
length of the path pe'- Then it is easy to see that the overall eost fix' of Too T' is 0{hk^/^) times that of tree 
T; this is due to the following implieations of the above lemma: (a) the total (buying) eost of all edges in T' 
is at most 2hk}/^ times that of T sinee eaeh edge is reused at most 2hk^^^ times, and (b) for any terminal 
X ^ X, the edges on its path to the root in T' eorrespond to disjoint sub-paths in the unique path between x 
and r in T, and henee the total length eost in T' is at most that in T. 

Using this lemma, we ean now eomplete the reduetion by “embedding” the tree T' in the layered graph 
Gh. 

Lemma 26. If there is a feasible solution of overall cost ffor the single-sink buy-at-bulk problem on G, then 
there is a feasible solution of overall cost O^hk^/^cff) for the same problem on the h-level layered extension 
Gh. 

Proof: Let T be the union of the paths in the optimum solution on the graph G. It’s easy to see that T 
is a direeted in-tree. First, we use Lemma 25 to transform T to tree T' of height h. As noted earlier, the 
overall eost fx' of T' is 0{hk^/^(l)). Now, we eonstruet a feasible tree Th in Gh using this solution T' as 
follows: eonsider eaeh edge {u,v) in T' where u is at level i and v is at level (i — 1). Then, inelude the edge 
{ui,Vi-i) in Th. Clearly, sinee T' eonneets all the terminals to the root, so does Th. Moreover, notiee that 
there is a 1-to-l mapping between edges in T' and edges in Th. 

To bound the objeetive value of the subtree, we relate the objeetive value for eaeh edge of the subtree in 
Gh to its eorresponding mapped edge in T'. First, note that the overall eontribution of an edge e' = {u', v') 
between layers i and i -\- \ towards 4>t' is equal to the sum of eosts and times the lengths of the 

edges on the assoeiated path pe' from u' to v'. This is beeause, by property (iii) of Lemma 25, the number 
of demands in the subtree rooted at u' is exaetly kf~^l^ and all of them traverse this edge to reaeh r. Next, 
we note that, by definition, the eost of the edge (tt', between layers i and i -|- 1 in T/j is equal to the 
shortest directed path from u' to v' in G aeeording to the metrie Cg -|- Sinee we ehose the shortest 

path, we get that the buying eost of edge G Th is at most the eontribution of {u',v') towards 

firp'. Moreover, the total length eost is at most times the length of the shortest path, whieh is at most 

the fixed eost of G Th (again, this uses the faet that there are exaetly terminals whieh 

route through this edge in Th also). It therefore follows that the overall eost of the solution that we have 
induetively eonstrueted in Gh is at most twiee the overall eost of T', whieh is at most 0{hk^^^)4>. ■ 

Theorem 7 follows from Lemmas 24 and 26. 
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